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ABSTRACT 



Aims. We re-examine the conditions required to steadily deviate an accretion flow from a circumstellar disc into a magnetospheric 
funnel flow onto a slow rotating young forming star. 

Methods. New analytical constraints on the formation of accretion funnels flows due to the presence of a dipolar stellar magnetic 
field disrupting the disc are derived. The Versatile Advection Code is used to confirm these constraints numerically. Axisymmetric 
MHD simulations are performed, where a stellar dipole field enters the resistive accretion disc, whose structure is self-consistently 
computed. 

Results. The analytical criterion derived allows to predict a priori the position of the truncation radius from a non perturbative 
accretion disc model. Accretion funnels are found to be robust features which occur below the co-rotation radius, where the stellar 
poloidal magnetic pressure becomes both at equipartition with the disc thermal pressure and is comparable to the disc poloidal ram 
pressure. We confirm the results of Romanova et al. 2002 and find accretion funnels for stellar dipole fields as low as 140 G in the low 
accretion rate limit of 10^'Mo.yr"'. With our present numerical setup with no disc magnetic field, we found no evidence of winds, 
neither disc driven nor X-winds, and the star is only spun up by its interaction with the disc. 

Conclusions. Weak dipole fields, similar in magnitude to those observed, lead to the development of accretion funnel flows in weakly 
accreting T Tauri stars. However, the higher accretion observed for most T Tauri stars (M ~ 10"*'Mo.yr"') requires either larger stellar 
field strength and/or different magnetic topologies to allow for magnetospheric accretion. 

Key words. Accretion, accretion discs - Magnetohydrodynamics (MHD) - Methods: numerical - Stars: pre-main sequence 



1. Introduction 

Classical T- Tauri stars (CTTS) are magnetically active, show 
evidence for circumstellar accretion discs, and can have mean 
photospheric ma gnetic field magnitudes around IkG (e.g. 
iJohns-Krull et aljri999l) . Such a strong stellar field is enough to 



disrupt the inner accretion disc, provided that one really mea- 
sures the large scale magnetic field and not only local strong 
multipolar components from starspots. However, this does not 
[ seem to be the case since recent polarimetric measurements 
(IValenti & Johns-Krulll 12004 ) indicate a weak dipolar compo- 
nent lower than 200 G. Moreover, observations show evidence 
for non-direct accretion. Inverse P-Cygni profiles with strong 
redshift absorption wings are indicative of polar accretion near 
free-fall veloc ities along magnetospheric field lines from the i n- 
ner disc edge jEdwards et alJll993HBouvier et al.lll999ll2003l) . 

Although t he magnetic field str ucture of these stars is prob- 
ably complex (Gregor v et al.ll2006l) . dynamical models of star- 
disc interaction usually assume an aligned dipole field, to sim- 
pUfy the analytical work. Under this assumption, stellar field 
lines threading the Keplerian disc below the co-rotation radius 
'"co - {GMt/Q.iy^^, where 0» is the stellar angular velocity, 
would lead to a spin-up of the star, versus those beyond r^o, 
to a spin-down. The radial extent of angular momentum ex- 
change between the star and the disc is then determined by (1) 
the disc truncation radius r^, where the magnetic dipole diverts 
the radially accreting flow to funnel flows; and (2) an outer ra- 
dius rout beyond which no more stellar field lines are connected 
to the disc. In this framework, a star-disc interaction occur- 



ring oi l a large radial ex tension (as p roposed by Ghosh & LambI 
19791: ICameron & Cam pbefl 19931 lArmitage & Clarkd 1 19961) . 
may lead to a disc-locking situation where the star remains at 
a slow rotation rate, despite accretion. On the other hand, it has 
been argued that this scenario is unlikely, since the stellar field 
lines would be ope ned up b y differential rotation until sever- 
ing this causal link (|AIv & K uiipers 1990; Lovelace et al.ll 19951; 
lMatt&Pudr"it3l2005 " for a recent discussion on that issue). The 
outcome of this latter scenario would be a star-disc interaction 
limited to a small radial extension around the disc truncation 
radius. Many theoretical models then assume that the disc in- 
ner edge should be close to the co-rotation radius (Konigl 19911 
IShu et alll 19941) for the sake of angular momentum equilibrium. 

For a decade, many nu merical works have investigated 
the star-disc inter action issue dHavashi et alJll996l ; [Kuker et akl 
120031: iLong et aT]|2005b . However, the formati on of accretion 
curtain s seems difficult to reproduce. Although iMiller & Stond 
(Il997h found such polar accretion in the case of a kG dipolar 
stellar field associ ated with a disc field in the same direction to 
the latter. .Romanova et a n (12002) were the first to demonstrate 
magnetospheric accretion along stellar field lines for kG pure 
dipolar field in axis ymmetric sim ulations and next by perform- 
ing 3D simulations dRomanova et al. 2003). 

In this paper, we address the issue of the disc truncation ra- 
dius and its localisation as a function of the disc (accretion rate) 
and stellar (dipole field) parameters. In Sect. 2 we provide ana- 
lytical constraints for driving steady accretion funnels and derive 
an estimate of the position of the truncation radius. We then use 
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numerical MHD simulations in Sect. 3 to verify this prediction. 
We confirm magnetospheric accretion for a slowly rotating star 
with an inner disc hole and a weak stellar magnetic field com- 
patible with observations of weak accretors. 



2. The disc truncation radius 

For a given accretion disc model, predicting where the truncation 
by the stellar magnetosphere will occur is an important issue. 
Different estimates were given in the literature. We consider here 
a pure dipolar field with a strength at the stellar surface in the 
equatorial plane equal to B,. 

A first dimensional estimate expresses the truncation radius 
Tt in the form rt - kr^ where the Alfven radius 



rA 



4pI2 



2GM,M„ 



1/7 



(1) 



is a characteristic length which can be derived by equating the 
ram pressure of a free-falling sp herical envelope wit h the mag- 
netic pressure of a dipolar field (Eisner & LambI 19771) . Different 
estimates of the adimensional coefficient k have been given in 
the li t erature, ranging f rom 0.5 (Ghosh & Lamb 1979; Konigl 



1991 



1995 



jLpng et al.l |2005D to ~ 1 (,Aronst,1993i: ,Ostriker & Shui 
Imng 1996). 

A second criterion (^Cameron & Campbeljl Il993t 
lArmitage & Clarkd 119961: iMatt & Pudritz 2005. and refer- 
ences therein) states that accretion funnels will take place 
when the magnetic torque due to the stellar field matches the 
"viscous" or turbulent torque, namely 
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where is the toroidal field at the disc surface, h the disc 

■ 1 

scale h eight, P the disc pressure and the lShakura & Svunvaevi 
( Il973h parameter This criterion (hereafter A) gives an upper 
limit for the truncation radius, as it only defines a radius where 
the star-disc interaction starts to affect accretion. This maximal 
truncation radius rt_n,ax occurs when the plasma beta, defined 
as j6 = 2fioP/Bl, becomes /3 - 2q/{ays) where q - IB^/B^I is a 
measure of the magnetic shear and s - h/r is the disc aspect 
ratio. In a thin Keplerian accretion disc, one gets /? » 1 at rt max 
sinc e q is close to unity to avoid opening of the magnetosphere 
(e.g. lMatt & Pudritzll2005h . 

A third criterion represents a more conservative approach, 
and gives only a lower limit rt,min. It states that accretion funnels 
take place when accretion is no longer possible because of the 
overwhelming field strength. It is usually written (iKoldoba et al.l 
120021: iRomanova et al.ll2002l) 



2 . 
— = pv^ -I- P 



(3) 



where v is the total speed. Note that this criterion is not predic- 
tive since we cannot calculate a priori the dominant azimuthal 
velocity v^ as it is itself an outcome of the star-disc interac- 
tion. Nevertheless, if one wishes to provide an estimate of the 
truncation radius, then we can use instead B^//io = pv^, where 



V0 = Vk = VGM»/r and a negligible thermal pressure (both ap- 
proximations are valid in a thin disc). This will be our criterion 
B. At the truncation radius estimated as such, the plasma beta is 
/3 = 2s^ « 1. 



From the previous discussion, it appears quite obvious that 
y6 ~ 1 is a better approximation for r^. The stellar magnetic field, 
which is bound to become dominant in the magnetosphere, must 
first favour accretion, i.e. the magnetic torque must be negative. 
If this is not the case, namely if r > > i"co, the disc mate- 
rial is radially expelled. This is the "propeller" regime as studied 
e.g. bv'Ustvu gova et all ( 12006 ) and references therein. Accretion 
thus implies rt < r^o, with stellar magnetic field lines as leading 
spirals. Below co-rotation, accretion will proceed quite naturally 
thanks to both the viscous and the stellar torques. There are then 
two more independent constraints that must be fulfilled to pro- 
duce steady funnel flows. First, the accretion flow must be pre- 
vented by the presence of the magnetosphere. The simplest way 
to express this is to require that the magnetic poloidal pressure 
balances the accretion ram pressure pu^. This defines a radius 
rb f where 



-2 



(4) 



where = Ur/Cs is the sonic Mach number measured at the 
disc midplane. Now, at radii r < rt^max, accretion is mainly due 
to the stellar torque and ms = 4q/yS » as: material is accreting 
to the star much faster than in the outer accretion disc. Second, 
material at the disc midplane must be lifted and loaded onto the 
stellar field lines. With a dipole field configuration, such a verti- 
cal motion can only be due to a vertical thermal plasma pressure 
gradient. It therefore requires that the magnetic field compres- 
sion is not too strong. This l eads naturally to an equipartition , 
yS ~ 1 (as already proposed bv lPringle & Reeslll972tlAlvlll980h . 

One important point to note is that once a large scale mag- 
netic field is close to equipartition in an accretion disc, it is 
able to deviate a large fraction of the disc plasma from its 
radial motion to a vertical one. This has been shown with 
the calculations of rnagnetized accretion-ejection structures by 
iFerreira & Pelletier (1995) and confirmed by numerical MHD 
simulations (Casse & Keppens 2002; Zanni et al. 2007). In one 
sense, making funnel flows involves the same physics as loading 
mass in magnetized jets. In fact, almost all of the disc mass can 
be lifted and loaded onto the field lines wh en ~ 1 , depending 
mostly on the field bending (see Fig. 3 in IFerreira & Pelletien 
1995). This is why we assume, for finding a simple analytical 
criterion, that the disc truncation radius r, is close to the radius 
where /3 ~ I, namely r, ~ rhf. Using the above estimates, we 
derive the following two constraints: 



/3 ~ I and nig ~ I. 



(5) 



that must be fulfilled in order to provide steady state funnel 
flows. For a dipole field, this translates into a theoretical trun- 
cation radius 



^ - 2 mfBlI'M-^I'M^'I'R'J' 

^ ci ^ ^ 



(6) 



where the stellar field B, has been normalised to 140 G, disc 
accretion rate Ma to 10~^ M© yr"', stellar mass to 0.8 Mq and 
stellar radius to 2R0. Note that these are all typical values for 
CTTS while the chosen value of the magnetic stellar field B, is 
consistent with observations of dipole fields in such objects (see 
Valenti & Johns-Krull 2004; Bouvier et al. 200 7, and references 
therein). It is clear that the conditions given by Eq. (|5]l state that 
the accretion speed close to the base of the accretion funnel is of 

the order of the Alfven speed. Indeed Eq. (|6]l can be rewritten as 

2/7 

'"t.th ~ '"A, where the Alfven radius for spherical accretion ta 
was defined in Eq. ([TJ- On the other hand, the co-rotation radius 



Bessolaz, N. et al.: Accretion funnels onto weakly magnetized young stars 



3 



is VcoIR, = 7.8 mI'^R-^pI'^ for a typical 8-day stellar period P„ 
which provides 



- 0.25 ml''Bt'^M-/'^M;'°'^'Rl^'^Pf" 

i * * ^ * 

' rn 



(7) 



Remarkably, taking typical values for CTTS with a low magnetic 
dipole field gives a theoretical truncation radius not only smaller 
than the co-rotatio n radius but a lso consistent with observations 
of inner disc holes (iNaiita et alj[2 003). This is in strong contrast 
with unobserved large scale kG fields usually used in the litera- 
ture. This implies that the formation of funnel flows should al- 
ways spin up the star unless (i) enough stellar field lines remain 
connected to the disc beyond rco and/or (ii) ejection of stellar an- 
gular momentum is indeed taking place. In the following section, 
2.5D numerical MHD simulations are used to show that the con- 
ditions in ^ are indeed those prevailing at the disc inner edge. 
If this is true, then we should recover a truncation radius located 
at our theoretical estimate (|6]l for a low dipole field of 140 G. 



3. Numerical experiments 

3.1. Equations and numerical setup 

We use the VAC Q code (lT6fl]|ll996l) to solve the full set of 
axisymmetric dimensionless resistive MHD equations (with the 
magnetic permeability yuo = 1) in cylindrical coordinates (r,z): 



^ -H V ■ (pvv - BB) -H V + P 



-pVOc 



dB 

dt 



H- V ■ (vB - Bv) = V X [77(V X B)] 



(8) 
(9) 
(10) 



dt 



£ + y -i-Plv-B(B-v) 



= 7772-B-(Vx77J),(ll) 



where p is the plasma density, Vp the poloidal velocity, v the 
total velocity, B the magnetic field, P the thermal pressure, rj 
the magnetic resistivity, J = V x B the current density, Oc - 



CM. 



the gravity potential created by the central star and E 



+ + ^ + P'^c the total energy density. The gravity is 
treated as a source term in the momentum and energy equations 
of VAC. 

Time evolution is done with a conservative second order ac- 
curate total variation diminishing Lax Friedrichs scheme with 
minmod limiters applied on primitive variables, except for den- 
sity where a van Leer limiter is used instead to better resolve 
contact discontinuities. Powell source terms are used to ensure 
the divergence free property of the magnetic field and the code 
has been modified so as to compute onl y the deviations from 
the dipolar component (Tanaka 1994; Po well et al.lll999l ). The 
splitting technique has three main advantages. First, it is crucial 
to properly represent an initial force-free configuration numer- 
ically. Second, since the total conserved energy contains only 
the energy associated with the deviation from the dipolar field, 
this method improves the computation of the thermal energy in 
low /3 regions. Third, when used in association with the Powell 
method, it helps control the divergence of the magnetic field. 



since only the deviation from the background field is used to 
calculate the divergence and the Powell source terms. The diver- 
gence method used here gives |V-B|A/? < 0.01|B| but simulations 
done with PLUTO0 and a constraint-transport scheme show that 
our results are not strongly modified by this method (Zanni et al, 
in prep). 

3.2. Boundary conditions 

Boundary conditions take symmetric and asymmetric conditions 
for the axis and the disc midplane, while continuous (outflow) 
conditions are used at the outer edge. At the inner edge cor- 
responding to the stellar surface, a linear extrapolation of den- 
sity and pressure is made. After using free extrapolation for the 
poloidal velocity, this latter is forced to be parallel to the to- 
tal poloidaljiiagiwticjield. The poloidal field is fixed, contrary 
to iRomanova et al] (|2002|) . We use a boundary condition on the 
toroidal magnetic field designed by Zanni et al. (in prep.). It al- 
lows us to derive by forcing the magnetic surfaces to locally 
rotate at the stellar velocity. In practice, the radial derivative of 
the toroidal field is computed from the angular momentum con- 
servation equation where the temporal derivative giving the local 
acceleration is replaced by 



p|(Q.) 



■p- 



(n,r + VpB^/Bp-ar) 
At 



(12) 



where the timescale At - AR/Va is the Alfven crossing time of 
one cell AR at the inner edge of the simulation box. This allows 
the system to evolve without imposing any arbitrary conditions 
on B^, and hence on the torques. 

3.3. Initial condition 

In the initial conditions we take a Keplerian disc surrounded by 
an adiabatic corona in hydrostatic (non rotating) equilibrium. 
The disc is adiabatic with an index y = 5/3 and an aspect ratio 
e - h/r - Cs/Vk = 0.1. The surface of the disc is determined by 
the pressure equilibrium between the disc and the corona, while 
the initial truncation radius has been chosen arbitrarily. The den- 
sity and pressure expressions for both the disc and the corona 
are 



Pd - PdO 
Pd 



2 

5? 



Rq 



-Il-.6^l^ 



f.1 +^2 



Pc 
Pc 



Ro 



PcO 
^^KOPcl) Pr 



(13) 



http://www.phys.uu.nl/~toth/ 



The density contrast between the disc and the corona is 
pco/Pdo - 0.01. A pure dipolar magnetic field is set up in the 
computational domain in equipartition with the thermal pres- 
sure of the disc at r = Rq. The rotation period of the star is 
set in order to place the corotation radius at r = 2/?o- The 
disc has no viscosity but is resistive with an alpha-like mag- 
netic resistivity decreasing on a disc scale height, namely v^, = 
Qfn,f2k/^^ex/?(-(j)"*). In this paper, we fix = 0.1. The initial 
poloidal flow within the disc is zero and the disc is slightly sub- 

Keplerian with = ^/l"^^fe2 
- A description of this code can be found in lMignone et alj (l2007h . 
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t= 0.0 



t= 5.1 



t=10.2 




t=15.3 



t=20.4 



t=25-5 




-3.^2 -.5.17 -2,8J -2,48 -2.14 -1.79 -1,45 -l,ig -D.76 -0.41 -0,07 



Fig.l. Resistive MHD simulation for a 5-day period CTTS with B, = 141G and = 0.1 after t=0,5. 1,10.2,15.3,20.4,25.5 
Keplerian periods at the disc inner edge corresponding to a physical time of 1.5 months. We show the density distribution in the 
computational domain using a log scale. The black lines draw the magnetic field lines and the black arrows represent the velocity 
field. The white line on the first snapshot represents an initial magnetic field line anchored at r^o- We superimpose a part of the 
computational grid to show the good resolution near the truncation radius. An accretion column is formed between rt and rtt (see 
text) and one observes the expansion of the poloidal magnetic field and transient disc ejecta. The accretion rate at the stellar surface 
is equal to 1.9 x lO^^Mo.yr"' at t=5 and stabilises around 0.91 x lO^'^Mo.yr"' at t=15. No X-winds are formed and the star is being 
spun up. 



The grid and stellar rotation period were chosen to allow 
good resolution at the truncation radius Rq while maintaining 
the co-rotation radius at r = 2Ro well inside the domain. Our 
polar grid of NrXNb - 170 x 100 is stretched in the spherical R 
direction (see Fig.[T]) and goes fromi?niin = 0.35/?o at the stellar 
surface to Rm^x - lQ35Ro- 



The results are presented in adimensional units: lengths are 
given in units of Rq, which corresponds to the truncation ra- 
dius of the reference simulation (see Sect. 3.5); speeds are ex- 
pressed in units of the Keplerian speed Vku = ^/GMTJRo and 
densities in units of pjo, which is the initial disc density at 
(r - Ro,z = 0). Time is given in units of the Keplerian pe- 
riod at 7?o, i.e. to = IttRq/Vko- Mass accretion rates is given in 
units of Mo = PoVkoRq while we express the torques in units of 
Lq = PoVl 



3.4. Reference values 

We consider a M* = O.SM© young star with a radius of Rt - 2Rq 
and a pure stellar dipole field with B, = 141G. With these 
assumptions the normalisation units will be Rq - 2.86/?» = 
4 X 10 V Vko = y/GMJRo ^ 1.63 x lO^m s"' and p^o = 
[BlRfj I {l^oe^Vl^Rfj - 5.51xlO"'°kg m"l Since we place the 
corotation radius at 2Ro, the rotation period of the star is f , = 



2n 



ylSRl/GM, 



5 . 1 days while the time will be scaled in units 

of to = 2nRo/VKO = 1 .54 x lO^s ^ 1 .78 days. The computational 
domain extends up to 0.3 AU. Finally the normalisation for the 
accretion rates is given by Mo = PoVkoRq - 2-27 x lO^^M© yr"' . 

From these reference values, simulations done here can be 
scaled for another range of parameters (/J*,Mt,B») in the follow- 
ing way: 
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'^U - "-"-^ \0.8Mo; \2RrJ ' 

Po = 5.51x10-"' kg 

5.1 (^f(^r'^^ days, 
.„.1.54xl0^(^f (^)-'^, 
Mo = 2.27 X 10- {^f {^f {^y"' Mo yr-. 



3.5. Methodology 

For these fixed stellar and accretion disc parameters, we make 
three simulations with different initial truncation radius r,,,. Our 
reference simulation (si) corresponds to an initial truncation ra- 
dius fulfilling Eq. namely r,o - 1 with y6(r,o) = 1 . Simulation 
(s2) is done for r,o = 1.5 with /3{r,o) = 10 and simulation (s3) for 
r,o = 0.5 with Pin,,) = 0.1. 

We then let the system evolve and observe whether or not the 
real disc truncation radius r, converges towards the theoretical 
radius r,//, as given by Eq. (|6]l. With our normalised quantities, 
the truncation radii in simulations (s2) and (s3) are thus expected 
to converge towards 1, with /3{rt) ~ 1. To check this, we identify 
this radius and then compute the plasma beta. It is however not 
straightforward to define this radius since the magnetic field is 
not a solid wall : the radial to vertical deviation of the flow is 
quite smooth. In practice, we get the truncation radius r, by de- 
tecting a steep decrease in density at the disc midplane (see Fig. 

3.6. Results 

Figure 1 shows a series of snapshots of our reference simula- 
tion (si). The upper left panel shows the initial condition. We 
have superimposed a part of the computational grid to show the 
resolution achieved. We have a resolution of 8 points in the ver- 
tical direction within the disc at each radius, while the resolu- 
tion within the accretion column reaches 20 points at the stel- 
lar surface. Magnetic field lines are in black while the white 
line traces a magnetic field line connecting the star to the co- 
rotation radius r^o = 2. After a rapid transient phase with the 
opening of stellar field lines, a quasi-steady situation is achieved 
where quasi-steady accretion columns are formed, even with a 
low stellar dipole field. The final (equilibrium) truncation radius 
is rt ~ 1.1 - 1.2, and thus well below the co-rotation radius. 

Figure |2] shows the projection F - Bp - F l\Bp\ of the various 
forces along a magnetic field line located at the centre of the ac- 
cretion column at f = 10. Since the stellar magnetic moment is 
directed northwards, a negative force is actually pulling material 
upwards. Accretion is achieved because of the negative magnetic 
torque (not shown here) but against the poloidal magnetic force 
Fm which tends to prevent it. Note however that it is negligi- 
ble with respect to the other forces. What drives the poloidal 
motion in the accretion column is actually the plasma pressure 
gradient Fp which was built in by the accumulation of accret- 
ing mass and allows it to be lifted up and loaded onto closed 
stellar field lines. It acts exactly as in accretion-ejection struc- 
tures, enabling the necessary tra nsition from the resistive M HD 
disc to the ideal MHD columns jFerreira & PelletiedlT995l) . The 
plasma pressure gradient remains dominant well above the disc 
surface, located around the curvilinear coordinate s - 0.2. Given 
the dipole topology, this is not surprising as material must first 
be lifted against gravity (Fc is initially positive) by Fp. Then, 



-0.40 t , , , \ , , , \ , , , \ , , , : 

0.0 0.2 0.4 0.6 0.8 



Fig. 2. Projection of the forces in normalised units along a mag- 
netic field line in the middle of the accretion column, for run (si) 
at f = 10. We represent the gravity Fq, the centrifugal force Fq, 
the thermal pressure gradient Fp, the poloidal magnetic force 
Fm and the total force F,,,! as a function of the curvilinear co- 
ordinate s. The disc midplane is located at 5 = and the disc 
surface corresponds to 5 ~ 0.2. Gravity begins to dominate the 
dynamics only at some distance between the disc and the star 
{s ~ 1.2 at the stellar surface). 

at some point (which depends mostly on the dipole geometry) 
gravity overcomes it and becomes the leading agent. Again, this 
is only possible because the centrifugal term Fc plays almost no 
role due to the azimuthal magnetic braking. Matter then reaches 
the star in a dynamical time scale with approximately free-fall 
velocities (vpoi ~ SlOkms"'). 

The strong differential rotation beyond rt, in a region where 
y6 » 1, leads to an expansion of the poloidal magnetic field lines. 
Such an expansion starts at the inner regions and enforces the 
outer field lines to inflate as well (but the cause there is not the 
differential rotation). Once these loops (inflated lines) reach the 
outer boundary of the computational domain, they open, mim- 
icking a reconnection. This opening is actually an effect of the 
boundary conditions used but, for all practical means, we see no 
strong bias on the evolution of the system. Anyway, on quite 
short time scales (f ~ 3), most of the stellar magnetic flux 
not related to the accretion funnels around n has been opened 
(Fig. 1). One might thi nk that a disc wind is drive n in this re- 
gion of the disc as the iBlandford & Pavnd (Il982l) criterion is 
fulfilled. Moreover, some mass is indeed leaving the disc along 
these open field lines. But this is only a breeze and not a proper 
jet: the ejected material does not reach super- Alfvenic speeds. 
We note also that no X-winds dShu et alJI 19941) are obtained de- 
spite the favourable magnetic configuration 0. The reason why 
no disc wind is obtained, neither extended nor X-wind, is that 
the field threading the disc in these regions is far below equipar- 
tition. 

The results of varying the initial truncation radius (simula- 
tions s2 and s3) are summarised in Fig. [3] In both cases, we 
observe a rapid convergence to rt ~ 1 . 1 - 1 .2 on a dynamical 
time scale. Although some fluctuations in time can be seen in rt, 

^ This is a direct simulation of the scenario invoked by these authors 
since the disc magnetic field comes from the opening of the magneto- 
sphere. 



6 



Bessolaz, N. et al.: Accretion funnels onto weakly magnetized young stars 




5 10 15 20 25 30 

Time 



Fig. 3. Top : Evolution in time of the position of the truncation radius r, (solid lines) for a set of resistive MHD simulations with 
Bt = 141G, Pt =5.1 days and - 0.1 and a different initial truncation radius. The time unit is the Keplerian period at r=l. All 
runs converge towards a truncation radius ~ 1.2 (solid lines). Dash-dot lines represent, for each simulation, the radius where 
j8 = 1 . We report also the radius rbf (dashed line) which indicates the base of the funnel flow, n.max (dashed line) and rt,min for our 
reference run (si). Bottom ; Evolution in time of the sonic Mach number at the base of the funnel, OTs(rbf), for run (si). 



it remains strikingly close to unity. Besides, it is really close to 
the radius of equipartition where yS = 1, ri, with a relative posi- 
tion lower than 4%. This is another indication of the importance 
of the plasma pressure gradient in defining the truncation radius 
and justifies the approximation y6(r,) ~ 1. 

The accretion rate onto the star is that which is actually 
observationally determined through, e.g. veiling measurements. 
It is obtained here by computing the mass flux in the accre- 
tion column, namely Ma - P^p ■ dS - -A-nRl pvR sin 9dQ. 
Surprisingly, it is found to converge towards 0.91 lO^^Mo.yr"' 
(see Fig. |5]l, hence a factor 10 smaller than the mean accretion 
rate in CTTS. We will come back to this issue later. If we now 
insert this value into Eq. |6] we find a theoretical truncation ra- 
dius rt_th = 1.38, using w., = 1. This is off by 20%, which is not 
bad considering our crude approximations of such a complicated 
problem. 

Since yS(r,) ~ 1 is well verified, the main source of discrep- 
ancy in Eq. |6] is due to the assumption of mj(r,) ~ 1 . This is 
too crude for an obvious reason. Indeed, the disc truncation ra- 
dius, as measured by the steep drop in density at the equatorial 
plane, is actually the point where Ur - 0, hence - 0. Figure[3] 
shows the evolution of the sonic Mach number m , computed for 
run (si) at rbf, which corresponds to the base of the funnel flow 
where the poloidal magnetic pressure matches the poloidal ram 
pressure. It can be seen that rt,f is larger than r, by almost 30%. 
This is not surprising as it is necessary to first brake efficiently 
material (rbt) before being able to lift it up {/,). This is also illus- 
trated in Fig. m where we plot the radial profile of several quan- 
tities at the disc midplane: density p, angular velocities (real Q 
and Keplerian Q^:) and nis- Clearly, taking nis = 1 to derive r, 
is too crude. We find that using a value of ^ 0.45, namely 
close to the real value msifhf) (see Fig.[3]l, provides a much bet- 
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Fig. 4. Radial distributions of density p, angular velocities Q. 
(real) and Q.k (Keplerian) and sonic Mach number nis - UrjCs 
at the disc midplane, for run (si) after 17 Keplerian rotations. 
Notice the large accretion velocity in the region of opened stel- 
lar field lines. 

ter estimate with rr,h close to r, with an accuracy better than 
10%. 

Although the simulations here do not take into account vis- 
cosity, for completeness, we plotted in Fig. [3] the evolution of 
the radius rt,max for run (si). To do so, we assumed - I (thus 
overestimating the viscous torque) and measured numerically q 
(see Sect. 2). As expected, rt,max remains always significantly 
larger than the real truncation radius rt, which implies that crite- 
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0.30 



Fig. 5. Evolution in time of the accretion rate Ma (dash dotted 
line), angular momentum flux transported by matter L„, (solid 
line) and by the magnetic field, computed for the closed Lf^ 
(dash line) and open Lfo (dotted lines) field lines. The units are 
normaUsed and the two hemispheres are taken into account. 



rion A is not good enough. On the other hand, criterion B would 
give a truncation radius located at r, - 0.3: no accretion col- 
umn should have been observed at all in our simulations. Our 
results clearly show that criterion B is not relevant. What about 
the criterion as expressed in Equation (O? We plotted in Fig. [3] 
the evolution of the radius rt^min given by this equation and com- 
puted for run (si). It turns out that it gives a good (though under- 
) e stimate of the real truncation radius r, as already pointed out 
bv lRomanov a et al. (2002). This is because the star disc interac- 
tion introduces a sharp decrease in both the disc midplane den- 
sity p and azimuthal velocity = Q.r (see Fig. HJl. However, 
as stressed in Sect. 2, while basically correct, such a criterion is 
meaningless as a predictive tool. 

As far as the forma tion of funnel flows is co ncerned, we fully 
confirm the results of iRomanova et all (l2002h : these flows are 
indeed robust features of axisymmetric MHD simulations. Their 
different boundary conditions on the magnetic field at the stellar 
surface, namely a fixed normal component and free conditions 
for Bg and B^, do not finally play any significant role for the 
truncation of the disc and the formation of accretion columns. 
We have also done simulations with other values of the mag- 
netic resistivity para meter a,„ with no change in the truncation 
radius. Note also that lRomanova et al.l (l2002h have included vis- 
cosity in their simulations while we did not, with no significant 
difference in the location of the truncation radius. To be more 
specific, using Eq. (|6]) with the higher accreti on ra t e of M g - 
9 X lO^^Mo.yr"' as measured by Romanova et al.l (l2002h . we 
find rt,th = 0.9 for - l.lkG using nis = 0.45, which is consis- 
tent w ith the truncation radius shown Fig. 16 in iRomanova et al.1 
(I2002h with a good accuracy. Fu r ther more, using Eq. ^ with the 
values provided bv lKuker et al.l (12003) one gets truncation radii 
smaller than the inner radial boundary, which explains why these 
authors did not find accretion columns. We are therefore confi- 
dent on our main conclusion, that is the validity of our criterion 



(©. 



field 
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dS = -AnRl pQ.r^VR sin ed6 



= r pQrhp-dS^-AnRl f 

Js Js 

= - I rB^Bp-dS=4nRl I rB^BR sin Ode 
Js Js 



(14) 



To explicitly write these expressions we assumed that the surface 
element is directed inwards with respect to the surface of the star. 
The flux carried by the magnetic field is the sum of that carried 
by closed field lines and open field lines, namely Lf = Lfc+Lfo. 
It actually corresponds to two possible electric circuits each re- 
lated to two different electromotive forces. A positive flux de- 
scribes a positive torque acting on the star and leading to a spin 
up while a negative flux leads to a spin down. The time evolution 
of these fluxes is shown in Fig. |5] for our reference run (si). Not 
surprisingly, the incoming angular momentum flux L„, due to the 
accreting material is positive but totally negligible with respect 
to that carried by the closed magnetic field lines Lfc, which is 
positive as well: the accreting star is only being spun up. When 
looking more closely at Lf, it turns out that there is a negative 
magnetic contribution to the torque due to the open field lines 
Lf„. This has also been previously reported in simulations and i s 
a natural outcome of the star-disc interaction dLong et al.ll2005h . 
We stress however that the actual torque is not controlled and one 
should not take it at face value. Indeed, this region of the mag- 
netosphere should be the locus of a stellar wind but the physics 
of its launching has not yet been addressed. Finally, let us turn 
back to our result that the accretion rate Ma measured onto the 
star is about 10 times smaller than expected. In fact, B, and M^ 
were considered in Eq. |6]as independent parameters. However, 
the mass inflow at the inner edge of the accretion disc is con- 
strained by the magnetic topology. Basically, the magnetosphere 
acts as a nozzle and the mass flow cannot be arbitrary. In a strict 
steady-state analysis it would be imposed by the regularity con- 
dition at the slow magnetosonic point. In our simulations, the 
flow reaches the slow point Msm = 1 above the disc at 5 ~ 0.3 
(see Fig. |6]l. Then, it remains supersonic with a maximum for 
the sonic Mach number around Mj ~ 3.5 but the funnel flow is 
always sub-Alfvenic with a maximum Alfvenic Mach number 
Ma ~ 0.3 near the middle path between the stellar surface and 
the disc inner edge. 



4. Conclusion 

We have confirmed that the formation of accretion funnel flows 
is a robust feature of axisymmetric star-disc interactions. We in- 
vestigated the physical conditions required to produce steady- 
state funnel flows onto a dipole and provided an analytic expres- 
sion of the disc truncation radius. We then used MHD simula- 
tions with VAC to show its validity. 

Our theoretical expression r, relates the disc truncation ra- 
dius to astrophysical parameters such as stellar dipole field B,, 
mass Mt, radius R, and accretion rate Mg . Although it resembles 
the Alfven radius ( lElsner & Larnbl 19771) sometimes invoked in 
(I " 



similar situations (^Bo uvier et al.l l2007h . it has been physically 



We now consider the fluxes of angular momentum, namely 
carried in by the infalling material and Lf by the magnetic 



motivated on very different ground. It is shown that it gives an 
accurate prediction of the real truncation radius as obtained in 
current 2D MHD simulations of a star-disc interaction. 

We report MHD simulations displaying accretion funnels 
with a weak stellar dipole field B, ~ 140 G from a resistive 
non-viscous disc. In this case, the disc inner edge is found closer 
to the star, below the co-rotation radius, in agreement with the 
size of inner disc holes. However, the accretion rate onto the star 
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Fig. 6. Sonic (dotted line), slow-magnetosonic (dashed line) and 
Alfvenic (dash-dotted line) Mach numbers within the funnel 
flow along the same magnetic field line at t=10 as in Fig.|2] 



which is imposed by the physics of magnetic accretion is mea- 
sured to be M ~ 10^"^ Mg.yr^K Even though th is magnitude of 
accretion rate is found for some T Tauri stars (see lGullbring et all 
Il998l) . this result shows that it is necessary to have stronger fields 
(by a factor of around 7) to have magnetospheric accretion for 
typical mean accretion rates of M ~ 10^^ M^.yr^^ . On the other 
hand, the few circular polarization measurements available, pro- 
vide an upper limit of 100-200 G for the dipolar component. 
However, recent spectro-polarimetric observations coupled with 
magnetic field reconstruction conducted on 2 CTTS V2129 Oph 
and BP Tau provide hi gher dipole field com ponents of 350 G 
and 1.2 kG respectively dPonati et al.ll2007aL Donati et al. 2007b 
in prep.). In any case, if the presence of kG dipole fields around 
CTTS are indeed ruled out, the requirement of forming quasi- 
steady accretion funnels imposes that the magnetic topology 
must be different. This is a firm result based on our dynami- 
cal calculations. One then needs to consider other stellar mag- 
netic field topologies such as multipolar components an d/or an 
inclin ed dipole configuration (Romanova et al. 2003; Lon g et al] 
|2007|) . In this paper, since there is no viscosity on our resistive 
accretion disc, we privileged a situation where rt to reach 

an open magnetic topology within the main part of the disc. This 
made it possible to have a braking torque feeding the disc even 
without viscosity. Apart from the localised magnetic flux chan- 
neling the accretion funnel, the other stellar magnetic field lines 
become causally disconnected from the disc. This forbids any 
disc-locking mechanism and the star is being spun up as a con- 
sequence of accretion. We also report the non development of 
X-winds despite the favourable topology. Both of these aspects 
deserve however a more detailed analysis as they may depend 
on the disc resistivity. While the disc truncation radius remained 
close to its predicted position for about 25 Keplerian periods at 
the disc inner edge, some variability is obviously taking place. 
This is a very promising topic as veiling measurements prob- 
ing ac cretion onto the star do show variability on different time 
scales encar & Basrill200d ). However making longer simula- 
tions requires one to implement the turbulent viscosity to allow 
for accretion within the disc beyond corotation. This work is in 
progress. 
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